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¢ Having read about “standard” PCKs and FKs in 
other tutorials you may want to learn about 
several “special” PCKs and FKs dealing with the 
Earth and the Moon. 


¢ While it is ultimately up to you, in most cases you 
should use the PCK and FK kernels described 
here when working with the Earth or the Moon. 


— If you need only “low accuracy” Earth or Moon orientation 
information, for instance for a visualization tool, then the 
standard text-style PCK data may be used. 
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¢ NAIF provides “high accuracy” orientation data for the Earth and 
the Moon in binary PCKs. 
— For the Earth, three versions are made: 
» High accuracy, frequently updated file 


* Contains high accuracy historical data and fairly accurate, short-term predict data 
» High accuracy, infrequently updated historical file 
» Lower accuracy long term predict file 


— For the Moon, a single, long-term file is made upon release of an official new JPL 
“Developmental Ephemeris” (DE). 


» Contains accurate historical and predictive lunar orientation data 


¢ To use these kernels: 
— Select binary PCK(s) having properties and time coverage that meet your needs 
» Unlike text PCKs, the time span covered by binary PCKs is limited 
Load the PCK(s) using FURNSH 
For the Moon, also load the lunar FK 
Reference the Earth body-fixed frame using the name ‘ITRF93’ 
» CAUTION: ‘IAU_EARTH’ cannot be used to reference high-accuracy Earth orientation data 


Reference a lunar body-fixed frame using one of these names: 
» ‘MOON_ME’ (Moon Mean Earth/Rotation axis frame) 
» ‘MOON_PA’ (Moon Principal Axes frame) 
» CAUTION: ‘I[AU_MOON’ cannot be used to reference high-accuracy lunar orientation data 
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Earth Binary PCKs 


High Accuracy Earth Rotation Model 
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¢ The ITRF93 high accuracy Earth rotation model takes into 
account: 
— Precession: 1976 IAU model due to Lieske. 
— Nutation: 1980 IAU model, with IERS corrections due to Herring et al. 
— True sidereal time using accurate values of TAI-UT1 
— Polar motion 


¢ It is more accurate than the IAU rotation models found in 
text PCKs. 


— See the plot on the next slide comparing orientation of the ITRF93 frame to that 
of the IAU_EARTH frame. 


» IAU_EARTH frame orientation error is ~0.06 degrees (~1 milliradian), or 
~6km on a great circle! 


¢ The highest accuracy is obtainable only for past epochs. 


— Unpredictable variations of UT1-TAI and polar motion limits the accuracy of 
predicted earth orientation. See plot on page 8. 
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Difference between the ‘I[AU_EARTH’ frame and the ‘ITRF93’ frame 


Prime meridian and pole offsets between ITRF93 and ISU_EARTH frames 
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Earth Predicted vs Reconstructed ITRF93 Plot 
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Difference between predicted and reconstructed orientation of ITRF93 frame 


Prime meridian and pole offsets between predicted and reconstructed ITRF93 frame orientation 
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¢ Binary Earth PCKs represent the orientation of an 
Earth ITRFxx body-fixed reference frame relative 
to the ICRF. 


— ITRF frames are defined by the International Earth Rotation 
Service (IERS). 


— Currently only the ITRF93 frame is supported within SPICE. 
» An update to a more modern version is planned. 


— Source data come from a JPL Earth Orientation Parameters 
(EOP) file. 


ICRF = International Celestial Reference Frame, often referred to in SPICE as the “J2000” 
frame, and also often referred to as the EME 2000 frame. This is an inertial frame. 


ITRF = International Terrestrial Reference Frame 
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¢ Three versions of the high accuracy binary Earth PCK are 
produced 

— “The latest,” using each new release by JPL of a reconstructed EOP file 

» Covers well into the past and approximately two months into the future 
beyond the production date 

» Accuracy of the future data degrades rapidly past the production date 
» Produced several times per week using an automated script 

— Long term predict, for uses beyond two months from the production date, and 

where the highest accuracy is not required 
» Produced infrequently 
» rhe several years into the past and approximately 30 years into the 
uture 


» Accuracy at epochs in the future is low compared to that for past epochs, 
but any of it is far better than what is obtained from the IAU rotation model 
for the earth (‘IAU_EARTH’) provided in any text PCK 


— History file, containing only high accuracy historical data 
¢ All are in the pck directory under generic_kernels on the NAIF 
server: https://naif.jpl.nasa.gov/pub/naif/generic_kernels/pck/ 
— Read the “aareadme’” file to see the file naming schema and more details 
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Accurate Earth Surface Locations 
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¢ High accuracy determination of surface locations 
relative to an inertial frame involves motions in 
addition to Earth rotation, including: 


— tectonic plate motion 
— tidal effects 
— relativistic effects 


¢ From the list of effects above... 


— Tectonic plate motion is accounted for in NAIF's DSN and some 
non-DSN station SPK files. 


» This helps ONLY station locations, not other surface 
objects. 


— The other two non-rotational effects affecting surface locations 
are not accounted for by any PCK, or by any other SPICE 
component. But the magnitude of these is under one meter. 


Special PCK and FK for Earth and Moon 


11 


Kernel Usage Summary: Earth 
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¢ To use high accuracy Earth orientation data 
— Load one or more binary Earth PCKs 


» If a long-term predict is used, load that kernel *before* 
loading any kernel containing reconstructed data so that 
lle pa data have precedence during the overlap 
period. 


— If your application uses any of the old, pre-N0062 APIs that 
make use of the default Earth body-fixed frame (see backup 
slides), load an Earth frame association kernel making ITRFxx 
the default Earth body-fixed frame. 


» But best to switch to use the “new” APIs that require you to 
specify which frame to use. 


* The new APIs are: ILUMIN, ILLUMF, ILLUMG, SINCPT, SUBPNT, 
SUBSLR 


¢ If you’re using SPICE to access Earth size and 
shape information, you'll also need to load a text 
PCK file containing these data. 
— Typically use the latest generic text PCK: pck000xx.tpc 
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Lunar Binary PCKs 


High Accuracy Lunar Rotation Model 
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¢ The high accuracy lunar rotation models available 
in binary PCKs are more accurate than the [AU 
rotation model found in a text PCK (‘TAU_MOON’). 
— For the time period of 2000-2020 the difference is approximately: 
» Worst case: ~0.0051 degrees, or ~155m on a great circle 
» Average: ~0.0024 degrees, or ~73m on a great circle 


— The error is due to truncation of the libration series in the [AU 
model 


— See the plot in the following chart comparing the IAU lunar 
rotation model to the integrated DE-421 model. 
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IAU_Moon vs MOON_ME Comparison Plot 
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Difference between the IAU_Moon frame and the Moon_ME frame (equivalent 
to the Moon_ME_DE421 frame) 


Prime meridian and pole offsets between MOON_ME and IAU_MOON frames 
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¢ Data for lunar orientation come from JPL’s 
DE/LExxx planet/lunar ephemeris files. 


— Binary lunar PCKs represent the orientation of the Moon’s 
“principal axis” reference frame-referred to as 
MOON_PA_Dexxx-relative to the ICRF*. 


ICRF = International Celestial Reference Frame, often referred to in SPICE as the “J2000” 
frame, and also often referred to as the EME 2000 frame. This is an inertial frame. 


JPL-produced planet/lunar ephemeris files are sometimes referred to as “DE/LExxx” but 
more often are referred to as simply “DExxx.” 
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¢ The high accuracy lunar orientation model obtained from the 
DE440 lunar ephemeris represents the result of a simultaneous 
numerical integration of lunar rotation and orbit, and of orbits of 
the planets. 
— The DE440 integration model includes*: 

» Relativistic equations of motion for Moon and planets. 

» Rotating Sun’s relativistic Lense-Thirring effect. 

» Solar radiation pressure on lunar orbit. 


» The gravity field of the Earth has zonal coefficients J, through J,;. The 
coefficient J, is a quadratic function of time. 


» Degree-2 tides on Earth depend on zonal, diurnal, and semidiurnal Love 
numbers with phase shifts for tidal dissipation. The diurnal and semidiurnal 
phase shifts depend on tidal period. 


» The degree-2 through degree-6 gravity field of the Moon comes from GRAIL. 
» The Moon has a solid lunar mantle and a fluid core. 


» oe. on the Moon are from the Earth, Sun, and planets Mercury through 
aturn. 


» Dissipation in the Moon depends on one tidal phase shift and one 
core/mantle interface parameter. The tidal Love number comes from GRAIL. 


» Torques on Moon from Earth J, interacting with lunar degree-2 figure. 
» Torque from an oblate lunar core/mantle interface. 


*Description provided by Dr. Ryan S. Park (JPL) 
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Lunar Rotation Model Effects-2 
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¢ The high accuracy lunar orientation model obtained from the 
DE440 lunar ephemeris represents the result of a simultaneous 
numerical integration of lunar rotation and orbit, and of orbits of 
the planets. 


— Parameters fit to lunar laser data for DE440 include’: 
» Initial conditions for the lunar orbit and orientation. 
» GM of Earth+Moon, which adjusts the mean lunar distance. 
» Lunar moment of inertia differences (C—A)/B and (B-A)/C. 


» Dissipation parameters for tides on Earth and Moon plus lunar core/mantle 
interface. 


» Locations of five laser retroreflector sites. 
» Thermal expansion and contraction of retroreflectors. 


— The best lunar range data of recent years have uncertainties < 1 cm. Our weighted 
rms residual is 9 mm. 


*Description provided by Dr. Ryan S. Park (JPL) 
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Lunar Frames Kernel 
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¢ A lunar frames kernel is available from NAIF. It has four 
functions. 
1. Make two lunar frames—Principal Axes (PA) and Mean Earth/rotation 
axis (ME)—known to the SPICE system. 
» Within SPICE their names are MOON_PA_DExxx and MOON_ME_DExxx 


» These frames are HN hig to a particular JPL-produced planetary and lunar 
ephemeris (e.g. DE421). 


» Special case for DE440: ME frame is MOON_ME_DE440_ME421, which is 
closely aligned with MOON_ME_DE421. 


2. Connect the MOON_PA_DExxx frame name to the high accuracy lunar 
orientation PCK data that implement the PA orientation. 


3. Provide specifications for implementing the rotation between the PA 
frame and the ME frame. 


» Makes the MOON_ME_Dexxx* frame available to SPICE. 


4. Provide generic frame names, aliased to the MOON_PA_DExxx and 
MOON_ME_Dexxx* frame names. 


» The generic frame names are simply MOON_PA and MOON_ME. 


» The generic names need not be changed in your programs when the 
MOON_PA_DExxx and MOON_ME_Dexxx* names change due to use of 
new defining data. 


» The DE-specific frames to which these aliases “point” will be updated 
by NAIF whenever a new binary lunar orientation PCK is produced. 
NAIF will release a new lunar FK at that time. 
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Kernel Usage Summary: Moon 
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¢ To use high accuracy Moon orientation data 
— Load the current binary lunar PCK 
— Load the current lunar FK 
— If your application uses any of the old, pre-N0062 APIs that 
make use of the default lunar body-fixed frame (see Backup), 


load a Moon frame association kernel making either MOON_ME 
or MOON _PA the default lunar body-fixed frame. 


» But best to switch to use the “new” APIs that require you to 
specify which frame to use. 
* The new APIs are ILLUMF, ILLUMG, ILUMIN, SINCPT, SUBPNT, SUBSLR 


¢ If you’re using SPICE to access Moon size and 
shape information, you'll also need to load a text 
PCK file containing these data. 
— Typically use the latest generic text PCK, such as pck00010.tpc 
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Lunar PCK/FK Summary 
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Which kernels are needed to access each of the three lunar 
body-fixed reference frames providing lunar orientation? 


Example of 
: oon pa de440 200625 
filename =» -P i = moon_220930.tf pck00010.tpc 

Kernel(s) : relate 

needed ex 
PCK 
Your 
objective 


Frame name to be 
used in SPICE 
software 


Usually a bad 
choice for the 
Moon! 22 
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Binary PCK File Format 


Binary PCK File Format 
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¢ SPICE binary PCK files are used to accommodate high 
accuracy rotation models. 


— Just as for SPKs and CKs, the data are held in SPICE Double Precision 
Array files (DAF) 
— Multiple types are supported 
» Type 2: Chebyshev polynomials are used to represent Euler 
angles giving orientation as a function of time. Rates are obtained 
by differentiating polynomials. Coverage intervals have fixed 
length. 
¢ Used for the Earth and the Moon 
» Type 3: Separate sets of Chebyshev polynomials are used to 
represent Euler angles and their rates. Coverage intervals have 
variable length. 
¢ Used only for Eros 
» Type 20: Chebyshev polynomials representing Euler angle rates 


¢ Provided to accurately represent “EPM” orientation data produced by the Russian 
Institute of Applied Astronomy (Ephemeris of the Planets and Moon) 


— Binary PCKs include a “comment area” for storing descriptive 
metadata 
» Access the comment area using the Toolkit’s commnt utility 
program 
— Binary PCKs support high-speed direct access 
» Chebyshev polynomials are fit to source Euler angles; these 


evaluate very quickly 
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Using Binary PCKs 


Precedence Rules 
for Text and Binary PCKs 
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¢ If two (or more) binary PCKs with functionally equivalent 
data are loaded, a later loaded file takes precedence. 


¢ Loading one text PCK that supersedes another can lead to 
errors if data from the “old” PCK remain in the kernel pool. 
— It’s essential to unload the old text PCK before loading the new one. 
» Use UNLOAD or KCLEAR to unload the old text PCK. 
— This problem doesn’t apply to binary PCKs. 


¢ If both a binary and a text PCK provide orientation for the 
same frame, data available from the binary PCK always take 
precedence over data available from the text PCK. 
— This is independent of file loading order 


— Note: the binary PCKs discussed in this tutorial define Earth-fixed and 
Moon-fixed frames different from those defined by a text PCK (e.g. 
pck00010.tpc), so there will be no conflict. 
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Tools for use with Binary PCKs 
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¢ Use the commnt utility to access a binary PCK comment 
area 
— Read, extract or insert metadata 


¢ Use the brief or spacit utility to summarize a binary PCK 
— brief is easier to use; spacit provides more information 


¢ Non-native binary PCKs can be read without first being 
converted to the native binary form 


— Converting a non-native binary PCK to native form will also speed up 
data access somewhat. Use toxfr/tobin or bingo. 


» toxfrand tobin are available in each Toolkit; bingo is available 
only from the NAIF website 
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Backup 


Association Frames Kernels 
for the Earth and the Moon 


Association FKs: Introduction 
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- In most SPICE modules that deal with one or more reference frames 
the name(s) of that/those frame(s) must be provided as input 
argument(s), for example: 

— CALL SPKEZR (target, time, frame, observer, correction, 
state, lighttime) 

¢ Many years ago the SPICE developers assumed there would be only 
one body-fixed reference frame associated with each natural body 
during a program run. 

— Thus a specific body-fixed frame name would rarely be needed as an input to 
modules dealing with body-fixed frames 

— Instead, SPICE could use the body-fixed frame associated with a given body simply 
by knowing the body name or ID 


» For most bodies SPICE associates the body with a body-fixed frame named 
IAU_<body name> (example: IAU_MOON) 


» This is known as the default body-fixed frame 
¢ This was a bad assumption... at least for the Earth and the Moon! 


— Multiple body-fixed frames exist for the Earth and the Moon 


— The default body-fixed frames for the Earth and the Moon, for which the defining 
data are provided in a generic text PCK (taken from an IAU report), are inaccurate 
(for the Moon) and very inaccurate (for the Earth) representations of the actual 
orientations of these bodies 
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¢ For the Earth and the Moon there are other 
choices for body-fixed frame that are almost 
certainly better than the default body-fixed frame 
conjured up by SPICE 


Body SPICE Default Body-fixed Frame Better choice 
Earth IAU_Earth ITRF93 


Moon IAU_Moon Moon_PA or 
Moon_ME 
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The Problem 
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¢ The SPICE modules that make use of the default body-fixed 
reference frame are these en Old: still available, but better to use those noted below 


— LSPCN, ET2LST, ILLUM, SRFXPT, SUBPT, SUBSOL (and their C, Icy and Mice 
equivalents) 


— Your code might overtly call one of these, or it could call one indirectly through use 
of a parameterized dynamic frame 
¢ NAIF rules regarding stability of our software offerings prevent us 
from changing the designs of those modules 
— So we have provided you means to change the default body-fixed frame associated 
with any solar system body of interest to you. See the next several pages. 
¢ However, starting with the version N62 Toolkits, a new set of 
modules is available for those calculations where precision body 
orientation is important. 
— These modules require the user to name the desired body-fixed frame, rather than 


using a default body-fixed frame ; 
New: safer to use, and offer improved 
— The new modules are these: pe accuracy in some cases 


» ILUMIN, ILLUMF, ILLUMG, SINCPT, SUBPNT, SUBSLR 
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Changing the Default Body-Fixed 
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¢ All bodies for which a body-fixed frame is defined by the IAU, and 
where the defining data are found in a SPICE text PCK file, have an 
associated default body-fixed frame name within SPICE: 
— The name pattern is: TAU <body name> 
— Examples: IAU MARS, IAU MOON, IAU_ EARTH 


¢ A different default body-fixed frame name can be assigned within a 
program by placing the following assignment in any text kernel 
that is loaded into the program: 

OBJECT <body name> FRAME = ‘<new default frame name>’ 
— Example: oBJECT MOON FRAME = ‘MOON ME’ 


¢ NAIF offers three “association FKs” to accomplish this. 
— See next page. 
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Using Association FKs to Change 
the Default Frame in the Older APIs 
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¢ For the Earth and the Moon, changing the default body-fixed frame 
name as described on the previous page can be accomplished by 
loading the appropriate “association” frame kernel provided by 
NAIF. The association kernels available are shown below 
— For the Earth: 
» Sartly BSso0C 1Urih93. tr 
— For the Moon: (pick one or the other—not both) 
» Moon assoc mest 
» Moon assoc pa.tf 
¢ These kernels are available on the NAIF server 
— For the Earth: 
» https://naif.jpl.nasa.gov/pub/naif/generic_kernels/fk/planets/ 
— For the Moon: 
» https://naif.jpl.nasa.gov/pub/naif/generic_kernels/fk/satellites/ 
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Lunar FK/PCK/Association FK Usage 
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Which additional kernel is needed to use the indicated frame in 
those (older) SPICE APIs* that use a default (assumed) frame? 
Pick one or the other. 


File name ==|_ moon assoc pa.tf moon_assoc me.tf 


PA ME 
Association Association additional 


FK** FK** kernels 


-LSPCN, £72157 EO SRERSUSERSUESOID 2) 


i < : replacements for these 
(and their C, Icy and Mice equivalents) four, which don’t use a 
default body-fixed frame: 
¢ ILUMIN, ILLUMF, ILLUMG 
**Any version of one or the other of these kernels is good indefinitely; : SUEONT 
you do not need to use the latest instance offered on the NAIF server. e SUBSLR 
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